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Abstract. We studied the evolution of two small flares (GOES class C2 and CI) 
that developed in the same active region with different morphological characteris- 
tics: one is extended and the other is compact. We analyzed the accuracy and the 
consistency of different algorithms implemented in RHESSI software to reconstruct 
the image of the emitting sources, for energies between 3 and 12 keV. We found 
that all tested algorithms give consistent results for the peak position while the 
other parameters can differ at most by a factor 2. Pixon and Forward-fit generally 
converge to similar results but Pixon is more reliable for reconstructing a complex 
source. We investigated the spectral characteristics of the two flares during their 
evolution in the 3-25 keV energy band. We found that a single thermal model 
of the photon spectrum is inadequate to fit the observations and we needed to add 
either a non-thermal model or a hot thermal one. The non-thermal and the double 
thermal fits are comparable. If we assume a non-thermal model, the non-thermal 
energy is always higher than the thermal one. Only during the very final decay phase 
a single thermal model fits fairly well the observed spectrum. 



1. Introduction 

Since February 2002, the RHESSI (Reuven Ramaty High-Energy Spec- 
troscopic Imager) satellite carries out X-ray and 7-ray observations of 
flares and other forms of solar magnetic activity (Lin et al., 2002). In 
particular, the 3-25 keV spectral band is of great relevance. It is com- 
prehensive of the thermal and non-thermal domain of the bremsstrahlung, 
giving us information about the energy deposition of the particles ac- 
celerated in the impulsive phase. The low-energy non-thermal electrons 
are especially crucial in small events and in microflares. Besides, a 
broadened emission line feature at around 6.7 keV, mainly due to Fe 
ions, and a weaker line feature at around 8 keV due to Fe and Ni ions are 
observable and can give information on electron distributions diff'erent 
from the ones obtained from the bremsstrahlung continuum. The X-ray 
band 3-25 keV is observed by RHESSI with unprecedented spectral 
resolution (1 keV) and imaging capability (2.3"). Moreover, when both 
attenuators are out, RHESSI has a much larger effective area than any 
earlier instrument, providing information on low- level energy releases 
with much higher sensitivity. 
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These RHESSI characteristics allow extensive studies of microflares 
related to both imaging and spectra (among others: Krucker et al., 
2002, Benz and Grigis, 2002, Krucker and Hudson, 2004, Hannah et al., 
2004, Liu et al., 2004). Images at low energy of these small events are 
generally reconstructed using the Clean algorithm available in RHESSI 
software (Hurford ct al., 2002, Schwartz et al., 2002). In fact, in this 
process of image reconstructions. Clean, together with Pixon (Metcalf 
et al., 1996) and Forward-fit (Aschwanden et al., 2002), converge to 
the best solutions for a simple source morphology and the photometric 
accuracy seems to be uniform in different energy ranges higher than 
10 keV (Aschwanden et al., 2004). It may be of interest to test the 
photometric accuracy and the consistency of the image reconstruction 
algorithms for more complex sources and lower energies. 

It seems quite established that at the peak time of a flare the ob- 
served spectra up to « 10 keV can be fitted with a single thermal model 
(T 10 MK) of the photon spectrum. At higher energies a good fit 
can be obtained adding a non-thermal bremsstrahlung emission with 
energy cutoff below 10 keV and a power-law index 7 between 5 and 
8, meaning that the spectra of small flares are usually softer than the 
ones of large flares. Benz and Grigis (2002) obtained a good fit of a mi- 
croflare spectrum at the maximum phase also adding a second thermal 
component with higher temperature (T~ 25 MK), but they reject this 
interpretation because this hotter component is no more present during 
the cooling and then the behaviour of the photons higher than 10 keV is 
more consistent with the standard flare scenario of precipitating non- 
thermal electrons. During the decay phase instead, the spectra can 
be fitted either with a thermal component only (Liu ct al., 2004) or 
with a non-thermal component with an index 7 ~ 10 (Krucker et al., 
2002, Hannah et al., 2004). Therefore it is very important to study the 
behaviour of the spectra of small flares during their evolution in order 
to investigate if the hardness of the spectrum changes with time, as 
can be observed in larger flares that show a direct correlation between 
the hard X-ray flux and the spectral hardness (Fletcher and Hudson, 
2002, Hudson and Farm'k, 2002, Grigis and Benz, 2004). 

In this paper we present a study on the morphological and spectral 
properties of two small flares (GOES class C2 and CI) in the energy 
range 3-25 keV. In particular we compare the image reconstruction 
algorithms to carefully check the characteristics of these two events at 
these low energies, when both attenuators are out, and we investigate 
how their X-ray spectrum evolves with time. 
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2. Observations 



An observing campaign coordinated between ground and space based 
instruments was planned to observe flare events, sampling the solar at- 
mosphere from the chromosphere to the corona. Monochromatic images 
at several wavelengths, in the continuum and within the Ha line have 
been acquired by means of the tunable Universal Bircfringcnt Filter 
(UBF) and the Zeiss filter at the Dunn Solar Tower of the National 
Solar Observatory (NSO) / Sacramento Peak. The field of view (FOV) 
is of about 150" x 150" with a spatial scale of 0.5" x 0.5" and a temporal 
cadence of few seconds. Spectra have been acquired with the Horizontal 
Spectrograph in 3 chromosphcric lines (Ca ii K, He i D3 and H7), with 
a temporal cadence of 4 s. Simultaneously, spectroheliograms of the 
same region were obtained in transition region and coronal lines with 
the Coronal Diagnostic Spectrometer (CDS) aboard SOHO (Harrison 
et al., 1995). During this campaign we observed the region NOAA AR 
10061 (N10W20) on 2002 August 11 from 14:00 to 19:00 UT. Two 
homologous flares (GOES class C2 and CI) developed in this region 
around 14:40 UT and 16:25 UT. In this paper we mainly concentrate 
on the diagnostic provided by RHESSI in conjunction with Ha data. 
The flare dynamics, using also spectra from NSO and CDS, will be 
analyzed in a later paper. 

The image of the whole active region at 14:41:47 UT (Figure 1), 
acquired with the HASTA telescope (San Juan, Argentina) in the blue 
wing of Ha line, shows the two ribbons developed during the first flare. 
Unfortunately the FOV of the UBF (indicated in Figure 1 by the white 
box) includes only the western ribbon of the flare. 

In Figure 2 we show the high resolution images obtained with UBF 
in the Ha line center during the two flares. Both flares have the 
characteristics of eruptive two-ribbon flares and are indeed triggered 
by the eruption of different portions of the same filament. However 
they exhibit very different morphologies. The upper panels show that 
the upper eastern part of the filament disappears at 14:14:57 UT. The 
emission maximum of the western ribbon is reached at 14:41:48 UT 
when two extended parallel ribbons developed in NE-SW direction (see 
Figure 1), at the two sides of the filament. MDI/SOHO (Scherrer ct al., 
1995) full disk magnetogram taken at 14:27:34 UT is used to determine 
the position of the apparent magnetic neutral line and the alignment, 
done using MDI and UBF continuum, is estimated to be within 1". The 
neutral line follows very well the filament and the ribbon positions hint 
at a system of large coronal loops (« 70 000 km) in East- West direction. 
Bottom panels show that at 16:00:16 UT the emission from the first 
episode has almost disappeared and the whole filament is visible. At 
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Figure 1. Image of the active region NOAA AR 10061 acquired by the HASTA 
telescope in the blue wing of Ha, approximately at the peak time of the first flare 
(West at right, North at top). Two bright ribbons are clearly visible. The white box 
(150" X 150"), that includes only one ribbon, indicates the FOV of the UBF/National 
Solar Observatory. The black boxes indicate the FOV used for RHESSI image 
reconstruction (see Sec. 3). 



16:31 UT very compact ribbons develop nearby the central part of 
the filament, suggesting the presence of stressed and lower loops. At 
17:13:31 UT the filament structure is again in the pre-flare conditions. 

Figure 3 shows the light curves of the two events observed in soft 
X-rays at 1 - 8 A by GOES, in the energy bins 3-6 and 12 - 25 keV by 
RHESSI and in Ha line center by UBF. We notice that GOES recorded 
a CI flare at about 15:50 UT which did not occur in the observed 
active region, and that GOES and UBF light curves are very similar 
for the two studied flaring episodes. RHESSI was in its night during 
the impulsive phase of the first event and acquired data only during its 
decay phase. For the second event RHESSI exited from its night just 
after the very first spike at 16:25:54 UT, well visible in the Ha curve, 
missing the beginning of the impulsive phase. The significant increase 
of the counts in the time interval 14:56 UT - 15:12 UT, at energies 
higher than 10 keV, is due to charged particle precipitation event and 
then the data in this time interval cannot be considered. RHESSI had 
both thin and thick attenuators out for the whole time interval. 
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Figure 2. Evolution of the two flares. Images have been acquired with UBF at the 
center of Ha lino. In the upper right panel the white line indicates the apparent 
magnetic neutral line obtained from MDI/SOHO data. At 14:41:48, the time of the 
maximum of the first flare, only one ribbon is within the FOV (see Figure 1). At 
16:31:09 UT, the peak time of the second flare, the ribbons follow the underlying 
filament and appear very close to each other. 

3. RHESSI Image Reconstruction 

In RHESSI software 5 different algorithms have been implemented to 
reconstruct images from the 'back projection' map, which contains 
sidelobes due to the modulation pattern of the collimators. Aschwanden 
et al. (2004) carried out a comparative study on the accuracy of the 
algorithms between 10 and 60 keV for a strong flare with two resolved 
sources. We did analyze the efficiency and the accuracy of different 
algorithms in the energy range 3-25 keV when both thin and thick 
attenuators are out and for two different source morphologies, i.e.: an 
extended and a compact flare. This energy range is crucial for the 
understanding of small flares and it might present different problems 
because of instrumental performances. 

3.1. Extended flare 

The minimum FOV to include all features of this complex-source flare 
is 128" X 128" (larger black box in Figure 1), centered on the point 
with solar coordinates (340,80), which seems to be the center of the 
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Figure 3. Light curves of the two flares observed by GOES (top), RHESSI (middle) 
and UBF (bottom). In the middle panel, the vertical dotted lines indicate the times 
at which both the RHESSI images and spectra have been studied. Notice the strong 
increase between 14:56 UT - 15:12 UT in the energy bin 12 - 25 keV that is due to 
a particle precipitation event (see text). 
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Figure 4- Extended flare: isocontours (20 %, 40 %, 60 %, 80 % of the peak intensity, 
indicated in the upper right corner) of the images reconstructed with four algorithms 
at 14:50:30 UT. 

X-ray emission. We reconstructed the images with a pixel size of 2"x 
2" considering a time interval At = 60 s which is short enough to give a 
low flux gradient at the considered time. We set off the time variability 
option (see later). The choice of collimators that contributed usefully 
to the definition of images has been accurate even if there is a certain 
arbitrariness. 

Collimator 2 cannot be employed for imaging below 20 keV due to a 
failure during the early phases of the mission (Smith et al., 2002). We 
made separate back-projection images for each other collimator and 
discarded those which show no obvious source. At the end, we decided 
to use collimators 5-8 (FWHM of the lower grid = 20.4"). Collimator 
9 (FWHM of 186") is excluded, as its low resolution does not add any 
new features, given the field of 128". Collimator 7 is also excluded in the 
3-6 kcV range because its energy threshold is 7 keV and its resolution 
is 3 keV. The reduced number of detectors used (only 3 over 9 in the 
3-6 keV range) could raise concerns on the results of this comparative 
analysis. However, we will show in the next Sections that the source 
parameters are not strongly affected by the number of detectors used 
in the reconstruction process. 

We compare the results obtained using four imaging algorithms 
available in the RHESSI data analysis software (Pixon, Forward-fit, 
MEM-Sato, Clean). We tried also the MEM- Vis algorithm but we could 
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Table I. Parameters obtained with different algorithms for the main source 
of the extended flare at 14:50:30 UT (see text for explanations). 



Algorithm 


Energy 






Ftot 






H/C 


Pixon 


3-6 


357.4,75.9 


10.1 


17362 


35.0 


20.9 


0.30 


Forw. Fit 


3-6 


359.0,73.0 


9.3 


18240 


38.2 


32.2 


0.19 


MEM Sato 


3-6 


358.3,73.6 


8.0 


16880 


36.0 


24.4 


0.38 


Clean 


3-6 


353.9,79.1 


4.3 


15000 


64.1 


46.8 


0.24 


Pixon 


6-12 


351.2,79.5 


0.45 


757.87 


36.0 


21.4 


0.16 


Forw. Fit 


6-12 


352.8,75.9 


0.43 


673.26 


37.4 


32.9 


0.05 


MEM Sato 


6-12 


356.7,73.3 


0.32 


735.80 


35.9 


29.2 


0.26 


Clean 


6-12 


355.5,74.0 


0.18 


678.92 


60.0 


42.9 


0.25 



not obtain the convergence for this method. Due to the complexity of 
the source, to obtain a good reconstructed image with Forward-fit we 
had to use 4 clhptical Gaussians. 

We considered the three energy ranges 3 - 6, 6 - 12 and 12 - 25 keV, 
but we had to exclude from our analysis the latter range, because the 
pileup effect is important for this flare (see Sect. 4.2). No correction for 
pileup is currently available for imaging, because the magnitude of the 
pileup varies on the same short timescales as the modulation, thus an 
exact correction for pileup would require an a priori knowledge of the 
image at all energies (Hurford, private communication). 

In Figure 4 we show the contours of the images obtained with the 
different algorithms in the 3-6 and 6-12 keV range at 14:50:30 UT 
(all RHESSI times refer to the center of the adopted integration inter- 
val). The other selected times provide very similar results and are not 
presented here. 

For each reconstructed image we give in Table I the parameters 
obtained in the 3-6 and 6-12 kcV energy ranges, that might help to 
assess the stability and the congruity of the different methods: 

— The peak position Xp and jp (arcsec), found by a parabolic fit of 

the nine pixels surrounding the image peak, and the peak flux on 
the detector Fp (photon s~^ cm~^ arcsec"^). 

— The total flux over the image on the detector Ftot (photon s~^ cm~^) . 

— The equivalent widths and Wy along the x and y direction give 
an indication of the linear dimension in arcsec and are defined 
following Aschwanden et al. (2004). 
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— The [Halo\/[Core] flux ratio H/C, where [Core] is the flux per 
unit surface in the part of the image inside of a circle centered on 
the peak emission with a radius ~ 1.5 times the resolution of the 
finest grid used, and [Halo] is the flux per unit surface in the part 
outside of it (Aschwanden et al., 2004). 

The peak positions agree within few seconds of arc, ~ 1-2 px, and 
the total fluxes agree within 20 %. The peak flux values and the linear 
dimensions Wx,y differ more than a factor of 2: Pixon and Forward-fit 
give similar results while Clean finds the lowest value for the peak flux 
and the largest for the linear dimensions. This is however an expected 
result, because Clean does not deconvolve the source from the Point 
Spread Function (PSF) of the instrument. However, the flux integrated 
over the source (F ~ Fp Wy) differs up to a factor of 1.5 and is not 
strictly conserved for different algorithms as in the case discussed by 
Aschwanden et al. (2004). 

The [Halo] /[Core] flux ratio indicates how much each method spreads 
the photons detected over the FOV outside the real sources. However, 
we stress that, especially in the case of a multiple-source flare, the 
[Halo] takes into account also photons coming by other real sources 
present in the FOV. Forward-fit, as implicit in the hypotheses used to 
reconstruct the image, obtains the lowest value for this ratio. MEM- 
Sato finds a source size comparable to Forward- fit and Pixon, but it 
obtains the highest [Halo] /[Core]. Thus, we can say that MEM-Sato 
leaves more photons out of the source. 

Finally, we investigate how the image reconstruction process is af- 
fected by the presence of unmodulated flux due to a combination of 
scattered flare photons within the instrument and a non-flare back- 
ground flux (Aschwanden et al., 2004). Pixon and Forward-fit have 
options to detect and remove this component. The corrected images 
obtained with Forward-fit are very similar to the uncorrected ones giv- 
ing variations of a few percent on the parameters reported on Table I. 
Instead, the parameters obtained with Pixon differ up to 40 % (if we 
exclude the total flux that is almost unchanged) giving very different 
results even in the peak position (Axp = —8", Ayp = -|-9") while the 
general morphology of the emitting region appears to be unchanged. 

The Clean algorithm has the possibility to correct the image for 
the residual map: this procedure strongly affects the total flux that 
decreases more than 20 % and then differs more than 40 % from 
values found with other algorithms. However, the [Halo]/[Core] ratio 
decreases as well as the flux integrated over the source. This means that 
the correction for the residual map subtracts photons from the field of 
view, but mainly from the halo. 



9 



Pixon 



Forward Fit 



MEM Sato 



Clean 




1 10 



o 



6 - 12 KeV 



310 340 370 



Figure 5. The same as in fig. 4 at 16:27:58 UT. 



3.2. Compact flare 

To reconstruct the images of this compact event we could use a FOV 
of 64"x64" (smaller black box in Figure 1), centered on the point 
with solar coordinates (340,80), and a pixel size of 2"x2". Since the 
light curve (Figure 3) shows rapid variations with two rise phases, two 
maxima and a decay phase, the time interval required for the image 
reconstruction should be short and we set At = 20 s. To create the 
images we used collimators 3-7 with a FWHM = 6.8"of the finest grid 
considered. Collimators 8-9 were excluded, because their FWHM (96" 
and 186" respectively) exceeds the FOV and collimator 7 was excluded 
in the range 3-6 keV, as in the previous case. This flare has a simpler 
structure and we were able to use 2 elliptical Gaussians to reconstruct 
images with Forward-fit algorithm. 

We examined the four imaging algorithms in different evolution 
phases of the flare with the time variability option switched off. The 
general behaviour of the algorithms does not change with the flare 
evolution and in Figure 5 we show as an example the images recon- 
structed in the first rising phase (16:27:48 UT - 16:28:08 UT) for the 
energy ranges 3-6 and 6-12 keV. In Table II the parameters obtained 
for these images are reported. 

Analyzing the parameters given in Table II we see that the peak po- 
sitions found with all algorithms agree very well within 2" , and the total 
fluxes agree within 10 %, but only if we exclude the Clean algorithm. 
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Table II. Parameters obtained with different algorithms for the main source 
of the compact flare at 16:27:58 UT. 



A 1 fT^Vl^" rl TVl 


Tr.noTTnr 

J^neigy 




rp 


F.. ^ 

■T tot 




Wy 




Pixori 


3-6 


344.4,72.6 


12.8 


2023 


10.1 


9.5 


0.11 


Forw. Fit 


3-6 


343.5,72.8 


15.0 


2121 


10.2 


9.9 


0.04 


MEM Sato 


3-6 


344.3,73.1 


13.3 


1885 


8.4 


7.1 


0.21 


Clean 


3-6 


344.0,72.7 


5.3 


724 


17.0 


11.5 




Pixon 


6-12 


347.7,73.4 


0.9 


77.21 


10.2 


9.7 


0.05 


Forw. Fit 


6-12 


346.5,72.3 


0.7 


67.37 


10.4 


9.7 


0.004 


MEM Sato 


6-12 


344.6,72.4 


0.8 


75.13 


9.7 


8.8 


0.08 


Clean 


6-12 


344.7,72.5 


0.3 


19.15 


15.7 


10.2 





Pixon and Forward-fit give very similar results for the peak flux and 
the linear dimensions while MEM-Sato and Clean give different values 
resulting in a flux integrated over the source a factor of 2 lower. The 
image obtained with Clean shows actually very high negative residuals 
out of the main source: consequently the total flux over the FOV is 
much lower than that obtained with other algorithms and moreover 
the [Halo]/[Core\ ratio cannot be reliably computed. Even if we use 
the available option in order to correct the image for the residual map 
wc cannot obtain more reliable results. As in the case of the extended 
flare the MEM-Sato algorithm gives the highest value of [Halo]/[Core\ 
ratio. Finally we use for Pixon and Forward-fit the options to remove 
the unmodulated flux: in this case the parameters change only of a few 
percent. 

For both the extended and the compact flare we did some other tests 
to verify the stability of the parameters found. Changing the pixel size 
from I" to 4" does not change the parameter values more than 1 %. 
A similar variation is obtained with the time variability option on and 
off, indicating that if the time interval considered is short enough to 
give a low flux gradient this option can be dropped. For the compact 
flare in the energy range 12 — 25 keV, when the counts are very low, 
we could not obtain the convergence of any algorithm if we consider 
a At = 20 s, but we should use At = 60 s. Within this interval the 
counts are highly variable and to reconstruct an image it is mandatory 
to set the time variability option on. 

We can conclude that in our two cases of small flares, the results ob- 
tained by different algorithms do not depend on the chosen energy band 
within the 3 12 keV range even if the detector 7 has been excluded 
only in the 3-6 keV band. This indicates that the reconstruction 
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algorithms are reliable even when only few collimators are useful to 
detect the emission source (sec Sect. 3.1). Different algorithms may 
find different problems in reconstructing the images, but the problems 
seem not to be related to the extension or compactness of the source 
but rather to the complexity of the region. 

MEM-Sato always gives the highest value for H/C as found by As- 
chwanden et al. (2004). However, in all the cases we examined it does 
not under-resolve the source, giving always linear sizes of the same order 
of the other algorithms. Pixon and Forward-fit give similar results but 
in the case of the extended flare the images reconstructed by Forward- 
fit do not seem to show the real complexity of the source. Thus we use 
Pixon to reconstruct the images for the analysis of flare evolution. 

4. Spectra evolution 

We analysed the spectral characteristics for the spatially integrated flux 
of our two flares as they evolve in time. The count spectrograms were 
created using the standard RHESSI software (updated April 2004). We 
used the front segments of the detectors in the 3 - 25 keV energy range, 
with the energy binning set to 1/3 keV, excluding detectors 2 and 7, 
which have lower energy resolution (Smith, 2002a; Smith, 2002b). The 
full spectrum response matrix has been used to calibrate the data, and 
the SPEX code (Schwartz, 1996; Smith et al., 2002) has been used 
for the spectral fltting. We considered energies above 4 keV, because 
the lower energy intervals are not well calibrated yet. We tried to fit 
the observed count spectrum with three different models of photon 
spectrum: 

— A single thermal bremsstrahlung emission whose best-fit parame- 
ters are the temperature T of the isothermal emitting plasma and 
its emission measure EM. 

— A thermal plus a non-thermal bremsstrahlung emission whose best- 
fit parameters are, besides T and EM of the isothermal emitting 
plasma, the power law index 7 and the normalization of the power- 
law at the normalization energy fixed at 10 keV. We assume a 
power-law index fixed to 1.5 for energy below the energy cut- 
off, that is determined as the one that gives the lowest value 
(see Sec. 4.3). We recall that is the reduced one, computed 
considering the numbers of free parameters. 

— A double thermal bremsstrahlung emission that assumes the pres- 
ence of two isothermal components and gives as best-fit parameters 
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Figure 6. (a) Count spectrum of the night background, (b) Count spectrum of 
the post-flare background. The spectra are significantly different only for energies 
£ < 8 keV. 



the temperatures T and T/j and the emission measures EM and 
EM^, for the cold and the hot component, respectively. 



4.1. Choosing the background 

RHESSI is a high-background instrument and the background selection 
in the spectra analysis is crucial in particular for small and low-energy 
flares. We considered two background spectra: the night background 
obtained averaging over few minutes before RHESSI exits the night for 
each flare, and the post-flare background determined between 16:55:30 
- 16:57:30 UT, well after the decay of any flare. The night background, 
which includes all the instrumental contributions, represents a sort of 
lower limit, while for the post-flare background we point out that the 
chosen interval is more than two hours (and, in RHESSI terms, more 
than an orbit) far from the flrst flare, and almost half an hour from the 
second one. 

In Figure 6 we compare the two background spectra. We can see that 
the post-flare background is at least one order of magnitude higher 
than the night background for energy lower than 8 kcV, while the 
spectra are almost identical at higher energies. The two peaks at around 
10.5 keV and 13.5 keV, present in both spectra, are probably due to 
the activation of the germanium detectors. We then expect the choice 
of the background will affect the flux spectra only at low energies. 
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Figure 7. Evolution of the extended flare during the decay phase. Images in the Ha 
line center with overplotted isocontours (40 %, 60 %, 80 % of the peak intensity) 
of RHESSI images in the 3-6 keV (top) and 6-12 keV {bottom), obtained using 
Pixon algorithm. 

4.2. Spectra of the extended flare 

In Figure 7 we show the evolution of the extended flare during its decay 
phase. The RHESSI images are obtained using Pixon algorithm with 
detectors 3-8 (excluding detector 7 in the 3-6 keV band, see Sect. 3.1 
) and Ai = 60 s. We could select all available collimators because 
Pixon decides which ones to keep and which to discard. The precision 
of the alignment between Ha and RHESSI images is estimated to be 
around 5". Ha images show the decreasing intensity of the ribbon. At 
14:50:30 UT the RHESSI emission maximum coincides with the loop 
footpoints outlined by the Ha west-ribbon and only a weak emission is 
present in the northern part of the FOV tracing the inferred large loop 
system connecting the two ribbons visible in Figure 1. Enlarging the 
RHESSI FOV we found that, from this time on, no X-ray emission is 
present in the eastern ribbon location. During the decay the intensity 
of Ha footpoints decreases as well as the peak flux of RHESSI images 
(variations are within a factor 10) and the X-ray maximum emission 
moves from the footpoints to the loop top. 

At least up to 14:53 UT the counts arc > 2000 count s~^ per detec- 
tor, which is the estimated threshold for the presence of pileup effect 
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Figure 8. 14.50.30 UT (a) Energy spectrum after the background subtraction with 
no pileup correction: the feature between « 10 — 15 keV is due to pileup. (b) The 
spectrum after the pileup correction. 



when both attenuators are out (Smith, 2002a). When pileup occurs 
multiple photons arc recorded as a single photon with an energy equal 
to the sum of energies of the individual photons. In Figure 8 we show 
the spectrum before and after the correction for the pileup. While 
the correction works properly below 12 keV, at higher energies the 
spectrum seems overcorrected. Usually the pileup is present for large 
flares and affects particularly the 20 - 50 keV energy range, and the 
error for its correction is estimated of about 20 % when the spectrum is 
fiat (Smith, private communication). In our energy range the spectrum 
is very steep and the error could be much higher. The clear presence 
of pileup in the spectrum prevents us to reconstruct reliable images in 
the 12 - 25 keV range. 

Hence we decided to fit the spectra in the energy range 4-10 keV 
using the same time intervals of the images. Contrary to expectations, 
the background choice does not strongly affect the resulting parameters 
even in the low energy part of the spectrum, probably because of the 
high counts recorded for this flare. In table III we give the parameters 
for the 3 models of the photon spectrum, obtained subtracting the post- 
flare background spectra (see Sec. 4.1). At 14:50:30 UT and 14:52:30 UT 
the single thermal model fitted the observed spectra giving a high 
of about 4 while the other two models fitted the spectra with low 



15 



Table III. Parameters obtained fitting the spectra with 3 models during the decay 
phase of the extended flare. EM, EM^, T and T^. are the emission measure and 
the temperature, respectively, of the cold and hot thermal components, 7 is the 
exponent of the power- law for the non-thermal component. The value is also 
given. 
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2.1 


15:14:00 


th. 


0.8 
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comparable values of about 2. For the non-thermal model we set 
the cutoff at 7 keV (see Sec. 4.3 for the discussion on the best choice 
of cutoff value) and we obtain a very steep power-law index 7 ~ 10, 
as found in the decay phase of some RHESSI microflares (Krucker et 
al., 2002, Hannah et al., 2004). For the double thermal model the hot 
component has a temperature of about 14 MK, compared with the 
9 MK of the 'cold' component. 

At 15:14 UT instead, when the footpoint emission is decreased and 
the contribution of the loop top increases, the single thermal model 
adequately fits the spectrum (T = 8 MK and EM = 0.8 x 10^^ cm"^) 
reaching a of about 1. This is coherent with the scenario of an erup- 
tive two-ribbon flare characterized by a global magnetic field disruption 
and processes of energy release long lasting after the impulsive phase. 
The magnetic reconnection process, responsible for these energy release 
episodes, is very fast in the beginning of the flare and slows down in 
later phases (Carmichael, 1964, Sturrock, 1966, Hirayama, 1974, Kopp 
and Pneuman, 1976). We can postulate that at 15:14 UT there are 
no more energy release episodes and thus, when the footpoints of the 
loop do not contribute significantly to the total emission, the plasma 
emission is only thermal. 

4.3. Spectra of the compact flare 

The RHESSI images are obtained using the Pixon algorithm with de- 
tectors 3-7 (the detector 7 is not included in the 3-6 keV band). In 
the range 3-12 keV we considered At = 20 s, while in 12 - 25 keV 
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Figure 9. Evolution of the compact flare. Images in the Ha line center with over- 
plotted isocontours (40 %, 60 %, 80 % of the peak intensity that changes within a 
factor 5) of RHESSI images in the 6-12 keV {top) and 12 - 25 kcV (bottom). The 
images are obtained using Pixon algorithm with At = 20 s and 60 s, respectively. 
At 16:32:53 no reliable image can be reconstructed in the 12 - 25 keV energy range 
due to the low counts. 



we had to consider At = 60 s (see Sec. 3.2). In Figure 9 we show the 

evolution of the compact flare beginning with the rise phase observed 
by RHESSI (times are indicated in Figure 3). The RHESSI contours are 
very similar in 3 - 6 keV and 6 — 12 keV energy bins and we show only 
6-12 keV and 12 - 25 keV contours. Ha images show that the bright 
features develop with time along the filament which is well visible before 
and after the flare (see Figure 2). Also the contours of RHESSI images 
in the 6-12 keV energy band change with time following the direction 
of the filament. In the 12 - 25 keV the peak emission moves along the 
filament and the contours, elongated in the North - South direction, 
probably include the loop footpoints and suggest the presence of low 
loops crossing the filament. 

We fit the spectra obtained with At = 20 s. Since no pileup effect 
is present in this flare, the data are reliable in the energy range 4-25 
keV. We fit the spectra changing the energy range from 4-10 keV to 
4-25 keV and we see that the found parameters change within 1 % 
while the decreases by a factor of two increasing the range of energy. 
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Figure 10. vs. the non-thermal energy cutoff value at different times during the 
flare. 

With the single thermal model the fit provides high values (~ 
6-8) with high residuals at energies > 10 keV. To fit the observed 
data adding a power-law non-thermal spectrum at first we analyze 
the output parameters changing the position of the energy cutoff from 
6 keV to 11 kcV. The found 7 index can vary up to 15 % with random 
residuals. The resulting seems to depend on the energy cutoff as 
shown in Figure 10 at different times during the flare. The minimum 
value is always reached with a cutoff between 7 keV and 8 keV, so we 
adopted 7 keV as the standard value for the cutoff. A posteriori, we can 
justify the use of the data in the 4-10 keV range to fit the spectrum 
and the setting of the cutoff of the non-thermal model at 7 keV also for 
the first flare, when the pileup effect makes unreliable the high energy 
data. 

In Figure 11 we show, as an example, the spectra at the time of the 
second peak (16:30:52 UT - 16:31:12 UT) fitted with different models. 
The broadened emission feature at around 6.7 keV, corresponding to a 
group of emission lines mainly due to Fe XXV and Fe XXVI, is also 
taken into account by the SPEX code in the fitting procedure. Phillips 
(2004) warned about the fact that SPEX computations are based on 
ionization and recombination data that overestimated the fraction of 
Fe XXV at T 10 MK. 

As expected, the choice of background for this very weak flare heav- 
ily affects the parameters of the cold thermal spectrum: the EM ob- 
tained subtracting the night background is approximately twice the EM 
obtained with the post-flare background while the temperature is on 
the average 0.5 MK lower. On the contrary the fitting parameters that 
characterize the power-law or the hot components are not substantially 
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Figure 11. Fitting of the observed spectra (crosses) at 16:31:02 UT with the 3 

considered models. The triangles indicate the background. The result of the fit 
is indicated by the solid line and the thermal component by the dotted line, (a) 
Single thermal model, (b) Thermal plus non-thermal model {dashed line), (c) Double 
thermal model (the hot component is indicated by the dot-dashed line). 

altered by the background as was also found by Qiu et al. (2004). 
In table IV we give the parameters obtained for the 3 models taking 
into account the post-flare background spectra (see Sec. 4.1). In the 
rise and peak phases of the flare the non-thermal spectrum has a 7 
of about 7 while in the decay phase 7 increases up to 10, as in the 
decay phase of the other flare. The high value of 7 reached when the 
emission decreases, indicates a softening of the emission with respect 
to the peak phase. If we assume a double thermal model, the hot com- 
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Table IV. Parameters obtained fitting the spectra with 3 models during different phases 
of the compact flare. EM, EM^, T and are the emission measure and the temperature, 
respectively, of the cold and hot thermal components, 7 is the exponent of the power-law 
for the non-thermal component. The value is also given. 
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ponent is always present during the flare and has a temperature T ~ 
20 MK (compared to 10 MK of the cold one) and an EM two orders of 
magnitude lower than that of the cold component. The temperature 
we found is anyway quite lower than the 35 MK indicated for the 
so called "superhot" component first identified by Lin et al. (1981) 
and afterward defined as the component giving rise to the Fe XXVI 
emission lines observed in large flares with BCS on Yohkoh (see e.g. 
Pike et al., 1996). 

4.4. Discussion 

An important result is that, for both flares, only towards the end of 
the decay phase the best fit can be obtained with a thermal photon 
spectrum while during the flare the best flt can be obtained, with 
comparable values, adding either a non-thermal component or a 
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hot thermal component. The two models are equivalent on this sta- 
tistical base and it is difficult to select the best model to interpret 
the observations. We try to use the Fe line feature at 6.7 kcV, well 
visible in our observed spectra, to choose between these two models. 
Phillips (2004) computed synthetic spectra using the CHIANTI code 
and found that this feature becomes visible above the continuum only 
when T > 10 MK (with an EM = 10^^ cm'^). Its equivalent width 
gives the best temperature diagnostic reaching a value of 3 keV for 
temperature T 20 MK. We measured the equivalent width of this 
line in the count spectra and we found a value always < 1.0 kcV that 
indicates a thermal plasma with T < 12 MK. We hence believe that 
the observed spectra could be better interpreted adding a non-thermal 
spectrum to a thermal one. 

From the fit parameters we can compute the thermal and the non- 
thermal contribution to the flare energy budget, assuming that the 
accelerated electrons lose their energy primarily through Coulomb col- 
lisions in the thick target approximation (Lin et al., 2001): 

Eth = SksTVEM V (la) 



Enon-th = 9.5 X 102^^(7 - l)B{j - 1/2,3/2)^^0 (lb) 

where ks is the Boltzmann constant, T and EM the temperature 
and the Emission Measure of the thermal part, V the volume of the 
source, B{a,b) the standard beta function, A the normalization factor 
of the photon spectrum, Eq the energy cutoff, and At the time interval. 
T, EM, 7 and A are derived from the fitting. We estimated V from 
the flare size (see Wx,y in Table II). For the extended flare we have 
observations only during the decay phase and we find Et^ 1 x 10^° erg 
and Enon-th ~ 1 X 10^^ erg. For the compact flare we find that both Et^ 
and Enon-th are approximately constant during the flare with Effi Ri 
5 X 10^^ erg and Enon-th ~ 5 x 10^*^ erg. Thus, the non-thermal energy 
is higher than the thermal energy by about a factor 10 not only in 
the rise and peak phases of the second flare but also during the decay 
phase of both flares. The high value of the 7 index in this phase, that 
indicates a softening of the emission, balances the low value of the 
normalization factor of the photon spectrum A. This means that the 
non-thermal electron energy is the energy supply for the heating of 
the thermal plasma. This result can thus be considered an a posteriori 
support to the non-thermal model choice. 
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5. Conclusions 



We use RHESSI data and ground based observations to study two 
small flares (GOES class CI and C2) developed in the same active 
region with different characteristics: the first is a two-ribbon extended 
flare implying high coronal loops while the other is a compact two 
ribbon implying low-lying loops. This is reflected in RHESSI imaging 
at different energies that outline the loop morphology in the two cases. 

We analyzed the accuracy of different algorithms implemented in 
the RHESSI software to reconstruct the image of the emitting sources, 
for energies (between 3 and 12 kcV) that arc particularly important 
when studying microflares. We found that all tested algorithms give 
similar results for the peak positions and the total flux on the considered 
FOV, as found by Aschwanden et al. (2004) for higher energies. The 
peak flux values and the linear sizes differ at most by a factor 2. Clean 
gives always the lowest peak and the largest size (does not deconvolve 
from the PSF of the instrument) while the other algorithms give very 
similar results. Clean can give additional problems even if the flare 
structure is simple. In fact for our compact flare it shows very high 
negative residuals out of the main source and then does not conserve 
the total flux over the FOV. The available option to correct for residual 
map does not help to solve the problem. MEM-Sato always gives the 
highest value for H/C, i.e. tends to spread more photons out of the 
source than the other methods, as found by Aschwanden et al. (2004). 
Pixon and Forward-fit do converge to similar results. Both algorithms 
have the option to detect and remove the unmodulated flux present 
in RHESSI data. Forward-fit finds a unmodulated component that 
slightly changes the source parameters for both flares, while in the 
case of extended complex flare, Pixon finds a component that strongly 
changes the source parameters. On the other hand, Forward-fit is not 
able to reproduce the real complexity of the source, hence we used 
Pixon (without background correction) to reconstruct the images of 
the two flares during their evolution. In fact, while the reliability of 
Pixon is not surprising, the most relevant result of our study is that 
each other method can meet significant problems, and special care is 
always required employing them. In any case, we outline the importance 
of comparing different algorithms when analyzing events in order to 
establish whether the reconstructed features are real or not. 

We studied the spectral characteristics of our flares by fitting the 
count spectra with 3 different models of photon spectrum: a single 
isothermal spectrum, an isothermal plus a non-thermal spectrum and a 
third one with two isothermal components. The photon spectrum seems 
to depend on the flare phase. The single thermal model is adequate to 
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fit the observations only during the very final decay phase. During the 
flare, instead, the single thermal model provides values a factor of 
10 higher than the one obtained with the other two models. The non- 
thermal and the double thermal fits are comparable, and the choice 
between them can not rely on statistics. The measured equivalent width 
of the Fe line feature at 6.7 keV is low compared to that expected for 
a plasma at 20 MK. Although further investigation is necessary, this 
may suggest that the hot thermal component is not relevant to the 
total emission. For the non-thermal spectrum we determined the energy 
cutoff value at 7 keV as the one that gives the lowest in the fitting. 
For the rise and peak phase we found a photon spectral index of 7 - 
8, quite high for standard flares, but not so uncommon in microflares 
(Krucker and Hudson, 2004, Hannah ct al., 2004), while during the 
decay phase 7 increases up to 10, indicating a softening of the emission 
with respect to the peak phase. With such a steep spectrum and a 
low energy cutoff the non-thermal energy contribution to the energy 
budget of the fiare is always higher than the thermal one and perhaps 
the non-thermal electrons provide the energy for heating the thermal 
plasma through collisional loss. 

As a "by-product" of our spectral study, we found that the pileup 
effect, usually observed during large flares, must be taken into account 
also for weak flares and that the correction enabled in the RHESSI 
software can give serious problems for the recovering of the spectra in 
the 3-25 keV energy range. We want to stress that the presence of 
pileup also affects the image reconstruction process, yielding unreliable 
images, because the photon energies are uncertain. 
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